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Abstract 

We present an algorithm for enumerating exactly the number of Hamiltonian chains on regular 
lattices in low dimensions. By definition, these are sets of k disjoint paths whose union visits each 
lattice vertex exactly once. The well-known Hamiltonian circuits and walks appear as the special 
cases k — and k — 1 respectively. In two dimensions, we enumerate chains on L x L square 
lattices up to L = 12, walks up to L — 17, and circuits up to L = 20. Some results for three 
dimensions are also given. Using our data we extract several quantities of physical interest. 

1 Introduction 

The subject of Hamiltonian circuits and walks plays an important role in mathematics and physics 
alike. Given a connected undirected graph G, a Hamiltonian circuit (or cycle) is a cycle (i.e., a closed 
loop) through G that visits each of the V vertices of G exactly once pQ. In particular, a Hamiltonian 
circuit has length V. Similarly, a Hamiltonian walk (or path) is an open non-empty path (i.e., with 
two distinct extremities) of length V — 1 that visits each vertex exactly once. Note that a Hamiltonian 
circuit can be turned into a Hamiltonian walk by removing any one of its edges, whereas a Hamiltonian 
walk can be extended into a Hamiltonian circuit only if its end points are adjacent in G. 

We add now to this list of well-known definitions the set of Hamiltonian chains of order k. 
Each member in Ck is a set of k disjoint paths whose union visits each vertex of G exactly once (see 
Fig. [I}. The set of Hamiltonian walks is then C\, and by convention we shall let Co denote the set of 
Hamiltonian circuits. Note that if V is even, Cy/2 is the set of dimer coverings of G. The Hamiltonian 
chain problem has been studied earlier by Duplantier and David on the Manhattan lattice [2], but 
never to our knowledge on an undirected lattice. 




Figure 1: Hamiltonian chain of order 4 on a square lattice of size 7x7. 



Figure 2: Transfer process for d — 2. The lattice is shown in solid line style, with the clotted lines 
representing lattice edges beyond the bounds of the lattice. The surface S is shown as a dashed line. 



Determining whether G contains a Hamiltonian circuit is a difficult (NP-complete) problem. An 
even more difficult problem is to determine how many distinct Hamiltonian circuits are contained in 
G. In this paper we shall present an algorithm that efficiently enumerates Hamiltonian circuits, walks, 
and chains for regular low-dimensional graphs. 

The motivation for studying such Hamiltonian structures is by no means limited to graph theory. 
Indeed, under appropriate solvent conditions, biopolymers such as proteins may fold to form compact 
conformations, the study of which is currently at the centre of an intense activity in the biophysics 
community. While real biopolymers contain complicated interactions which can probably not be fully 
accounted for within any simple lattice model, the study of Hamiltonian walks has been advocated as 
a first approximation for understanding qualitatively the excluded- volume mechanisms at work behind 
such problems as polymer melting [3] and protein folding [I]. Our extension to Hamiltonian chains 
permits to study polydisperse models of several polymers. 

Another interest stems from the study of magnetic systems with O(n) symmetry in physics. These 
can modelled on the lattice as self-avoiding loops (each having the weight n) [5], which, in the limit 
of vanishing temperature T, are constrained to visit all the vertices [6]. Coupling such systems to a 
magnetic field H amounts, in a perturbative expansion around H = 0, to inserting pairs of loop end 
points 7 . In the limit n — > 0, the partition function Z of the 0(n) model at T = in a weak magnetic 
field H can thus be expressed in terms of the number of Hamiltonian chains as 



Finally, the exact enumeration of configurations is useful for settling issues of ergodicity when 
developing algorithms that provide unbiased sampling of Hamiltonian walks in two [8] and three [9] 
dimensions. 

In section [2] we present our enumeration algorithm and discuss some aspects of its implementation. 
Results in dimensions d = 2 and d — 3 are given in section [3l For convenience, we limit the discussion 
to the simplest lattices (square and cubic) , although the construction extends straightforwardly to any 
regular lattice. Our results are strongest in d = 2 where we determine all Ck for L x L square lattices 
up to L = 12, C\ up to L = 17, and Cq up to L — 20. In d = 3 the largest lattice that we were able to 
access has size 3x4x4. Finally, we show in section |4] how to extract physically interesting quantities 
from our data. 



We first present our algorithm in dimension d — 2 and then discuss the necessary modifications for 
d = 3. For convenience, we limit the presentation to the simplest lattices, viz. an L\ x L2 square lattice 
and an L\ x L2 x L3 cubic lattice, although the construction extends straightforwardly to any regular 
lattice. The boundary conditions are free (non-periodic), although it will be clear that it is easy to 
introduce periodic boundary conditions along one of the lattice directions. 

The algorithm is based on the transfer matrix principle, according to which the lattice is cut into 
two parts by means of a conveniently chosen d — 1 dimensional oriented surface S. The part of the 
lattice above (resp. below) S is called the future (resp. the past). The surface S cuts the lattice only 
at mid points of edges. 




(1.1) 



fc 



2 Algorithm 



2 



Figure 3: Configuration on a partially constructed 6x6 square. 



At the initial (resp. final) step of the enumeration the whole lattice belongs to the future (resp. 
past), so the algorithm consists in sweeping S over the entire lattice. This is done by gradually pushing 
S towards the future, so that in any one step of the algorithm a single vertex is transferred from the 
future to the past. A few subsequent steps for a 4 x 4 lattice are shown in Fig. [2j 

In any step, the configuration of the system is described by some information about the edges cut 
by S, and by the number k of chains which have already been completed. The information refered to 
is the connectivity of the cut edges with respect to the part of the lattice which belongs to the past, 
and is best illustrated by an example (see Fig. [3} . The complete description of the configuration reads 
in this case (0, 1, 1, 2, 3, 2|0|I), where the first L\ entries refer to the state of the cut edges which are 
parallel to the 2-direction (vertical), and the next entry refers to the state of the one cut edge which 
is parallel to the 1-direction (horizontal). The last entry is the number k = 1 of completed chains. In 
the connectivity part of the information, we use the following coding: 

1. A zero entry means an empty edge. 

2. Two equal positive entries mean a pair of edges which are connected in the past by part of a 
chain. Each of these edges will eventually have to be linked to a chain end point in the future, 
so as to form a complete chain. 

3. An unpaired positive entry means an edge which forms part of a partially completed chain, one 
end point of which has already been fixed in the past. The edge will eventually be linked to 
another end point in the future, leading to the formation of a complete chain. 

It is important to avoid any redundancy in this connectivity information. A unique coding is obtained 
by requiring that the positive entries (i.e., unpaired entries, or the leftmost member of a pair of equal 
entries) be arranged in increasing order (1,2,3,...) when reading through the coding from left to right. 

In step t of the enumeration, the configurations are transfered from "time" t to time t + 1 . More 
precisely, each configuration at time t is examined in turn, and all its descendent configurations at 
time t + 1 are generated by exhausting the possible arrangements of the chain at the vertex which is 
transfered from the future to the past. The information described above is necessary and sufficient for 
deducing the connectivity information at time t + 1 from that at time t. Note that since the transfered 
vertex is not allowed to be empty, there are four possible arrangements if it accommodates a chain 
end, and six arrangements if a chain passes though it. 

Each configuration generated at time t + 1 is inserted in an appropriate date structure — a hash 
table — which also keeps track of its weight (here an integer). The weight of a descendent configuration 
is the sum of the weights of all the parent configurations that generated it. Some of the generated 
descendent configurations are however rejected before insertion in the hash table (see below). Once 
step t has been completed, the hash table storing the configurations at time t is erased, and one can 
move on to step t + 1. In this way, only two hash tables (at times t and t + 1) are needed in the entire 
process. 

Note that the choice of data structure is essential for the feasibility and the efficiency of the al- 
gorithm. A hash table permits to store the configurations via a key which is obtained by reading its 
coding as one large integer, modulo a suitably chosen prime. Storing and retrieving configurations can 
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be done in constant time, i.e., independently of the number of configurations being stored in the hash 
table. 

The hash table also allows to keep track of the weight of each configuration, according to the above 
rule. Namely, when a descendent configuration is generated with weight w, we first make an attempt 
of looking it up in the hash table at time t + 1. If it is not there, it is inserted with weight w. If it is 
already there, w is added to the weight of configuration already present. 

If in the transfer process the two ends of the same chain (coded by two equal positive entries) join 
up, the resulting configuration is rejected, since this would mean forming a cycle rather than a chain. 
(We make an exception to this rule when the very last vertex is added, since this permits to enumerate 
Co-) If an unpaired positive entry gets left behind in the past it means that a chain has been completed, 
and so k — > k + 1. 

Denote now a general configuration as (s2,i, S2,2, • • • , S2,Li|si|fc). At step t — 0, the initial state 
is (0, 0, ... , 1 0) and has weight 1. When a row of the lattice is completed, any configuration with 
s\ ^ gets rejected. When transfering the i'th vertex of the last row, any configuration with S2,i ^= 
gets rejected. This trick allows us to avoid having to deal with a lot of special cases when a boundary 
vertex is transfered — and also makes it much easier in practice to implement the algorithm correctly. 
After step t = LiL 2 the lattice belongs completely to the past, and all the configurations are of the 
form (0, 0, ... , 0|0|fc). Their respective weights are precisely the Ck that we wanted to compute. 

The maximum lattice size that we can attain is essentially limited by the number of different 
intermediate configurations generated in the transfer process. This number attains its largest value 
after transfering the next last vertex in the next last row. In practice we could store at most ~ 10 8 
configurations. 

As usual in enumeration studies, the coefficients Ck are much larger than the integers which are 
usually represented by a computer (< 2 32 for an unsigned integer on a 32-bit machine). We therefore 
repeat the enumeration several times, computing each time the result modulo different coprime integers 
(2 32 , 2 31 — 1, 2 31 — 3,. . . ), and reconstruct the true result in the end by using the Chinese remainder 
theorem. Note that the use of modular arithmetics is possible because the weights of configurations 
are constructed only by successive additions of positive integers. 

The counts for systems of size L\ x L 2 and i 2 x L\ should of course coincide. Verifying that this 
is indeed the case is however a very strong check of the algorithm, since permuting L\ and L 2 (with 
L\ 7^ L-i) leads to a completely different transfer process in terms of the propagation of the surface S. 
We have performed such checks both for d = 2 and d = 3. 

We now describe briefly how the algorithm can be adapted to a d-dimensional hypercubic lattice of 
size L\ x L 2 x • • • x Ld- The surface S is pushed though the lattice by means of d nested loops, of which 
the innermost loop (at nesting level d — 1) moves S along the 1-direction, etc., and the outermost loop 
(at nesting level 0) moves S along the d-direction. In general the loop at nesting level I moves S along 
the (d — ^-direction. 

A configuration is given by ({sd}|{s<i-i}| • • • |{ s i}|fc), where the space {se} describes the edges cut 
by S which are parallel to the ^-direction and consists of JIi=i entries. When the loop at nesting 
level d — I is executed for the last time, the entry in {se} corresponding to the position of the loops at 
nesting levels > d — i must be zero; otherwise the configuration is rejected. 

At each vertex there are 2d possible local arrangements if the vertex contains a chain end, and ( 2 2 d ) 
arrangements if it does not. 

We have implemented the algorithm for d — 2 and d = 3. It is of course most efficient in low dimen- 
sions when the number of entries necessary to describe a configuration is small, and the configurations 
themselves are strongly constrained by topology. We were however able to obtain useful results as well 
for d — 3 (see below). 

3 Results 

3.1 Two dimensions 

We first present our results in two dimensions. 
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_fc Ck_ 

1076226888605605706 

1 3452664855804347354220 

2 4105040900990127258563352 

3 1716401559105599779087093260 

4 363652056217217171206243035340 

5 46148041957435926988244999692732 

6 3870392655399966034741749180958852 

7 229000444797839686805213214595470648 

8 10014026193777241299766692880686035774 

9 335126632781634776435981605808153310160 

10 8818298873873444121262995871243309826506 

1 1 186425336415902384343216389461927330172268 

12 3222564357088784934867009058887596660853042 

13 46216822292126998413476281985324396507245748 

14 556693787783862608927984698386363470187981938 

15 5690759125611797657956588062161464969526622268 

16 49812855339352875263449851394183743992177894504 

17 376259523557799790490076563092079957367971441020 

18 2469068810121023544317004188530710352090945411914 

19 14159258603854781892361610528528465090872582482860 

20 71328606843660293099526723817890226149324996051770 

21 317097933449802440642304292096463208253752002853220 

22 1249084471204154623059161399602853755323161718519772 

23 4375441585007318378166364769735565649135514407725428 

24 13673227330703731913694049643952596405766113534103766 

25 38227041578177782993926722638883248316303238609485072 

26 95854543129195970883751499759209798433629811140560618 



27 216054324315261909273142199235822975041539861478121508 

28 438604247185012333557530043782600945755207568230643594 

29 803320261 178156739715403163647291125863198470203538048 

30 1329413609161899137884978107235777594261127805779481138 

31 1990427307505969944020236847124806666399811554129750496 

32 2699118726122396972411170136925985393885382039591724116 

33 3318040996252249925481740087801704155429556174343699608 

34 3700359327341155923554786919768922871490843218081792948 

35 3745875480760985174397430937750746156033853657449987776 

Table 1: Number Ck of Hamiltonian chains of order k (i.e., consisting of precisely k chains) on a 12 x 12 
square lattice. 
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k 



36 3443370502558807290562354000836385732268208485521139030 

37 2875008789367922659848083522243722363355224468273444220 

38 2180493763475645200773617970678658343878870756101785154 

39 1502095391370671575772987130076864642955725212507397512 

40 939647724819716693388579329116456350135042650270872340 

41 533562853839154855341927965097210146851707481332264436 

42 274862642496699290252290990160279057175930352324091824 

43 128362555399743719434877352579592630895860764794774624 

44 54294968510865162278386828132399087185256262942414932 

45 20778050812329434163278524326813392470106157931284728 



46 7184786845784650425312709323260198674084909532776292 

47 2241482341086337065894490074920307755810306693310936 

48 629823815715953133331560488991904014678562091268042 

49 159077245534404479393324552931624000438466735966020 

50 36035273248074550220433523348607010222299409542106 

51 7302543915028881257299361934945060798679283726808 

52 1320074215106942124873659331692519827058055405644 

53 212172409136637900106629576229565529242562898136 

54 30209940418516563152074763182692186029320843454 

55 3794668249506100568268876660836383536523653816 

56 418512747239130986282686885840526291044063692 

57 40310331424536568202368529058344353430213252 

58 3369925505474814797946214848424889586336952 

59 242794327043218279242923744053680409745316 

60 14951909347540490900851120307089218715790 

61 77949364740771636380822990073456271 1356 

62 34012688096545799769214778699372499928 

63 12253433732360441 19139709333819722548 

64 35847082988635315875798770248538788 

65 834223151185015104742566527439696 

66 15044012072281443196476714201368 

67 203123129951305502546186136496 

68 1958918610759335996516705296 

69 12602334728369293472453184 

70 48539905585658564517760 

71 91904378228899701504 

72 53060477521960000 



Tabic 1: (continued). 
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L Cj_ 

2 4 

3 20 

4 276 

5 4324 

6 229348 

7 13535280 

8 3023313284 

9 745416341496 

10 730044829512632 

11 786671485270308848 

12 3452664855804347354220 

13 16652005717670534681315580 

14 331809088406733654427925292528 

15 7263611367960266490262600117251524 

16 662634717384979793238814101377988786884 

17 66428994739159469969440119579736807612665540 

Tabic 2: Number C\ of Hamiltonian walks on an L x L square lattice, up to size L = 17. 

_± Co_ 

2 1 

4 6 

6 1072 

8 4638576 

10 467260456608 

12 1076226888605605706 

14 56126499620491437281263608 

16 65882516522625836326159786165530572 

18 1733926377888966183927790794055670829347983946 



20 1020460427390768793543026965678152831571073052662428097106 

Table 3: Number C\ of Hamiltonian circuits on an L x L square lattice, up to size L = 20. 

We have been able to solve the full Hamiltonian chain problem for L x L square lattices up to 
size L — 12. The number of circuits Co vanishes when L is odd by an easy parity argument, but the 
remaining Ck with k — 1,2,..., [L 2 /2J are all non-zero. The complete result for L — 12 is shown in 
Tabled] 

When L is even, C^i / 2 should be the number T>l of dimer coverings of an L x L square lattice. We 
have checked that our data agree with the analytical results [TU] for T>l for all L = 2, 4, 6, 8, 10, 12. 

If only the first few Ck are needed, the enumeration can be taken to larger sizes by rejecting all 
states (s2,i,S2,2, ■ • ■ ,S2,Li\si\k) with k > fc max . 

In particular, we have obtained the number C\ of Hamiltonian walks up to size L = 17; see Tabled 
This extends the L < 7 results by Mayer et al [Tl] by ten new terms. Note also that Jaeckel et al [12 
have proposed a Monte Carlo method for estimating C\ for larger L. These authors obtain 1.3582 • 10 7 
for L = 7 — that is 0.3 % above the exact result — and 2.7791 • 10 9 for L — 8 — that is 29 % below the 
exact result. 

Variant Hamiltonian walks, constrained to have their end points on diametrically opposite corners 
of an L x L square with L even, have been studied in [13] . Since this is technically an easier problem 
than our unconstrained walks, the enumerations could be taken to size L = 34. 

Finally, we have obtained the number Cq of Hamiltonian circuits up to size L = 20; see Table |3] 
This extends the L < 16 data [14] by two new terms. 
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k 2x2x2 2x2x3 2x3x3 



3x3x3 



3x3x4 




1 
2 
3 
4 
■5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 



G 
72 
204 
108 

9 



22 
584 
4204 
7604 
4541 
852 
32 



324 

16880 2480304 

270756 104844792 

1281376 1246834176 

2507084 6520250088 

2281064 17852434656 

985354 27873228036 

190580 25864316448 

13834 14445212196 

229 4798350687 



911288760 
91325100 
4119048 



57048 



7803954743412 
41553510978656 
134709106959932 
280608712776492 
388267754276278 
363680422635635 
232420898624633 
101132591631452 
29594655770318 
5683316575620 



3918744 
677849536 
40656040968 
816646740296 



687432832414 
48991382300 
1837669320 
29304199 
117805 



Table 4: Number Ck of Hamiltonian chains of order k on various cubic lattices of size L± x L2 x L3. 
Blank entries are zero. 

3.2 Three dimensions 

Our results for the full Hamiltonian chain problem on small L\ x L2 x L3 parallelepipeds are given 
in Table 0] In addition we find for the 3x4x4 system a number of Co = 3777388236 circuits and 
Ci = 1073054619800 walks. 

Note that the transfer matrix method is essentially limited by the area of the smallest cross section 
of the parallelepiped. It would thus be possible to extend the enumerations to some systems with, say, 
L\ = L2 = 3 and L3 > 5, but we have chosen to focus here on closc-to-cubic shapes which are the most 
challenging. 

The number of walks C\ has been much studied in the area of protein research [T51, \W\ , whereas the 
numbers Ck with k ^ 1 have to our knowledge not been considered previously. Note that the works 
[15l [16] were based on direct enumeration, meaning that in contrast to our transfer matrix method 
each individual conformation was actually generated. The limitation of direct enumeration is thus the 
number of conformations being counted, since the CPU time requirement is (at best) proportional to 
this. Accordingly, Ref. [TB] uses massively parallel supercomputer facilities to access the 3x4x4 
system. By contrast, our transfer matrix approach is limited rather by the memory than the CPU 
time. Unfortunately, this memory limitation put the 4x4x4 system just a little outside our reach. 
On the other hand, the counts for the 3x4x4 system were made very fast, in just a few minutes. 

In Refs. p~5j [16] the counts were produced modulo the symmetry group of the lattice. For the 
2x2x2 and the 3x3x3 systems, our results for C\ come out as exactly 24 times those of [TS]. For 
the 3x3x4 and 3x4x4 systems, our results for C\ are precisely 8 times those of [16 . This means 
that for each of these systems, all Hamiltonian walks are unrelated by lattice symmetries. 

4 Applications 

The enumerations reported above conceal many quantities of physical relevance. We discuss here some 
of them. 
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4.1 Conformational exponents 

The radius of gyration R of a polymer of length / 3> 1 is expected to scale like 

R~l" (4.1) 

where v is a standard critical exponent [7]. For Hamiltonian walks on an L d hypercube in d dimensions, 
we have obviously R ~ L and / ~ L d , and so 2/ = 1/d. Non-trivial information is however contained 
in the number of circuits and walks, here both supposed to have one marked monomer attached to a 
fixed point: 

Cx-VP" 1 , e ~n l r vd . (4.2) 

Here fj, is the so-called connective constant, and 7 is another critical exponent |T. In our setup, the 
circuits are unmarked and the end points of the walks are free to be anywhere on the lattice, and so 

C1/C0 - r +1 - L (7+1)d . (4.3) 

In addition to this leading behaviour there are subdominant corrections due to surface effects. 

In two dimensions, we can extract results for \i and 7 using the data in Tables GH3] Due to parity 
effects, this is best done by working in terms of the ratios ^nx^ < where r)(L) — Cq(L) for even L or 
Ci(L) for any parity in the case of /i, and rj(L) — C\{L)/Cq{L) for even L in the case of 7. The naive 
approximants are further extrapolated using standard finite-size scaling techniques. This gives 

H = 1.473 ±0.001, 

7 = 1.042 ±0.003. (4.4) 

Our estimate for \i is in good agreement with (but less accurate than) the currently best known 
estimate [B] 

l-i = 1.472801 ± 0.00001 . (4.5) 

Note that the latter uses exact predictions from field theory for the leading finite-size corrections, a 
scheme that we have not adopted here. The constrained Hamiltonian walks considered in [13j led 
to a consistent value for /j,. The long-standing history of numerical and analytical estimates for fi of 
Hamiltonian walks can be found in the introductions of [T2J. IS] . 

Our estimate for 7 is a nice confirmation of the exact field theoretical result [6] 

7= — = 1.04464 ••■ (4.6) 
'112 v ; 

Previous numerical results, as discussed in [6] and references therein, were obtained in a cylindric 
geometry and assumed certain results of conformal field theory The present estimate therefore furnishes 
a more direct verification of the exact result (14.61). 



4.2 Contact probabilities 

As already mentioned in the introduction, a Hamiltonian circuit can be turned into a Hamiltonian walk 
by removing any one of its edges, whereas a Hamiltonian walk can be extended into a Hamiltonian 
circuit only if its endpoints are adjacent. This implies that the probability that the two end points of 
a walk are adjacent is 

Padj = ■ (4-7) 

Li 

In two dimensions, the p a dj for L x L square can be computed from Tables HH3] The resulting 
numerical values are displayed in Tabic [5] Note that p a dj vanishes for odd L (as it does on any 
bipartite lattice having an odd number of vertices) . 

The values in Table [3] have been used in to test that a certain Monte Carlo algorithm for 
producing Hamiltonian walks did indeed give unbiased results. 
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8 


0.09819321919798768694 


10 


0.06400435120127304050 


12 


0.04488610346836087660 


14 


0.03315398616246303584 


16 


0.02545282308230371747 



Table 5: Probability p a dj that the two end points of a Hamiltonian walk on an L x L square lattice are 
adjacent. 



Physically, on may argue that p a( jj is proportional to the probability that the two ends of an open 
polymer (walk) join so as to form a ring polymer (circuit). 

It is tempting to try similarly to construct from the ratio of C\ and C2 the probability that the 
conformation of two chains is such that one end point of each are adjacent on the lattice. Unfortunately, 
this is not possible, since certain two-chains (more precisely those in which an end point of one chain 
is adjacent to both end points of the other chain) can be obtained from more than one one-chain by 
removing an internal edge in the latter. 



4.3 Lee- Yang zeroes 

The study of phase transitions through the location of partition function zeroes in the complex magnetic 
field plane was initiated by Yang and Lee [17]. In particular, these authors established that for the 
Ising model these zeroes lie on the unit circle in terms of the variable x = e 2H , i.e., they correspond 
to purely imaginary values of the field H. The zero closest to the positive real axis is denoted e 1,6 ", 
where 8 C is the so-called Lee- Yang edge. At the critical temperature, 9 C — > in the thermodynamic 
limit, and its finite-size scaling permits to access a critical exponent. Another possible approach is to 
study the density of zeroes g(0). Creswick and Kim [18] have shown that at the critical temperature 
g (0) „ \0\i-/t for <C 1. 

We have studied the zeroes of the partition function (|1.1[) in terms of the variable y — —H 2 for 
L x L squares with L < 12. For odd L there is one trivial zero at y = (since Co = 0), and for even 
L the two zeroes closest to y = form a complex conjugate pair with an imaginary part that tends to 
zero as L — > 00. Disregarding these "exceptional" zeroes, all the remaining zeroes (for L of any parity) 
are found to lie on the positive real axis in the complex y-plane, corresponding to purely imaginary H 
as in the Lee- Yang theorem. 

We now define a finite- L approximation to the density of zeroes in the point y n as g(y n ) = y + l- y ; 
where we have arranged the zeroes of Z(L x L) in increasing order j/i < J/a < . . . < yN- 

The approximations g{y) are shown in Fig. [4] for L = 9, 10, 11, 12. One observes a clear crossover 
near y = 1, separating two regimes of power law behaviours. For even L the curves bifurcate for y -c 1, 
which can be remedied by regrouping the zeroes two by two (not shown). The power laws extracted 
from the largest available sizes read 



.-0.46 



for y <g; 1 



We have no satisfying explanation for these exponents at present. The naive application of the 
standard scaling laws vd = 2 — a (Josephson) and a + 2(3 + 7 = 2 (Rushbrooke), and the results v = | 
and 7 = for the critical y — system [6], leads to l/<5 = — ^fg = —0.0218 • • • , which is clearly off 
the mark. 
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Figure 4: Density g(y) of partition function zeroes in the variable y — —H 2 for the Hamiltonian chain 
problem on L x L squares. 
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